Metastable Congested States in Multisegment Traffic Cellular Automaton 
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We investigate a simple multisegment cellular automaton model of traffic flow. With the introduc- 
tion of segment-dependent acceleration probability, metastable congested states in the intermediate 
density region emerge, and the initial-state dependence of the flow is observed. The essential feature 
of three-phased structure empirically found in real-world traffic flow is reproduced without elaborate 
assumptions. 
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I. INTRODUCTION 

The cellular automaton model has established itself as 
a standard in the analysis of traffic flow with its sim- 
plicity, affinity to computer simulations, and also with 
its extendibility to handle the realistic traffic 0, 0] • A 
particularly appealing aspect is its capability to capture 
the essential characteristics of the traffic flow with very 
few system parameters, which is a surprising economy of 
assumptions. With further assumptions for the adoption 
to the specific cases, the model is shown to be capable 
of reproducing down to the fine details of traffic flow in 
various settings, and it has been used for the analysis of 
real-life traffic with great success in past decade Q. 

Lost in this process of elaboration, however, is the el- 
egance of generality which the original standard model 
possesses. There is an obvious generic feature of the traf- 
fic flow that is beyond the reach of original model, but 
appears to await the explanation with a simple model. 
That is the three phase structure. In a typical traffic, the 
low and high density regions show the free motion and 
traffic jam. In the fundamental diagram of traffic, that 
plots the flux against the density, they are respectively 
represented by the increasing and decreasing lines, both 
of which are reproduced well with the standard model. 
In the intermediate density, however, empirical data in- 
dicate the existence of the third phase, in which there 
is no clear single line in the fundamental diagram, but 
the points are scattered in broad region, indicating the 
existence of the fluctuation in the traffic flow. 

In this article, we intend to show that this third phase 
appears in a very simple extension of the standard model 
that includes segments which has different system param- 
eters. These segments are introduced to represent the 
natural bottlenecks - narrow road segments, hilly seg- 
ments, a segment under construction, etc. - that are 
ubiquitous in real-life traffic. We base our discussion pri- 
marily on numerical simulations and supplement them 
with simple intuitive physical arguments. 

It has been recognized for some time that the fluc- 
tuation in the traffic flow is closely related to the the 
appearance of the metastable states 0, IE @- We find 
that the emergence of the third phase is also a direct re- 



sult of the metastable states that appear with a minimal 
assumption in our model. 



II. TRAFFIC MODEL WITH PROBABILISTIC 
ACCELERATION IN MULTI-SEGMENTED 
ROAD 

We first lay out our model. A road is divided into a 
one-dimensional array of L cells on which M cars are 
placed. The location (i.e. cell number) and the velocity 
of the cars are represented by integers Xi and Vi with the 
car index i — 1,2,..., M. Cars are indexed in ascend- 
ing order to the direction of the motion without over- 
lap; x\ < X2 < •■■ < xm- The road is considered to 
be circular, so that the (L + l)-th cell is identified with 
the fist cell. The road is subdivided into S segments 
Si = {x\0 < x < L x }, S 2 = {x\Lx < x < L 1 +L 2 }, 
Ss = {x\Li + ... + Ls-i < x < L}, with the constraint 
on the sub-length L\ + ... + L$ = L. Each segment, 
s = 1,...,iS, is assigned its own maximum velocity U s 
and non-acceleration probability R s . The time evolution 
of the system in discrete time steps t — > t + 1 is described 
by the following three consecutive updating for Vi and Xi 

1) Probabilistic acceleration : 

Vi — > min(fi + 1, U s ) with probability(l — R s ) 

2) Breaking to avoid collision : 

Vi — > min(wi, x i+ \ - x t - 1). (1) 

3) Advancement : 



The updating is done in parallel for all cars i — 1,2, ...M : 
Namely, the current position of the preceding car is used 
to calculate the distance Xi+\ — Xi — 1. 

The model can be thought of as a variant of the stan- 
dard Nagel-Schreckenberg (NS) model in that it follows 
the acceleration and deceleration step-by-step. In fact, 
if we change the ordering of "breaking" and "random 
deceleration" in the standard NS model, we obtain a 
model where cars accelerate, and then randomly deceler- 
ate. The model then becomes almost identical to ours, 
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apart from the treatment of the cars running at the max- 
imum velocities. This inverting of the order effectively 
prohibits the cars running in the jam from random decel- 
eration, thus models the "heightened alertness" of drivers 
who have the preceding cars in their close sight. Our 
model is also related to the Fukui-Ishibashi (FI) model 
in that both employ probabilistic acceleration, which is 
in contrast to the NS model that treats the process by 
splitting it into acceleration and probabilistic decelera- 
tion. 

Using the second step of the rule (QJ , we may define a 
jam as a block of cars in which the distance of neighboring 
cars is below the maximum velocity, A crucial feature 
of our model is the existence of the metastable block of 
jam made up of equi-spaced cars. This is a direct result 
of our model in which no further random deceleration 
occurs for cars inside the jam. The only limiting factor 
of the growth of this block is the percolation at the head 
of the block whose rate is proportional to (1 — R s ). This 
feature captures oft observed self-forming bottleneck of 
slow-moving cars in single-lane roads in real-life traffic. 
This could be regarded as a main justification for devising 
yet another variant of traffic cellular automaton model. 



III. NUMERICAL RESULTS 

We perform numerical simulations on the model pre- 
sented above with one segment, two-segments, and six 
segments cases, the last one being the model of generic 
realistic traffic. 

The macroscopic characteristics of traffic systems are 
best summarized with so-called fundamental diagram, 
which is the traffic flux F plotted as the function of traffic 
density p, each of which are defined as 



i=l i=l 

Here, () t stands for the average over many time steps t. 

We first show the results of the single-segment traffic 
in FIG. 1. The size of the road is set to be L = 200 
and the maximum velocity is chosen to be U = 8. Fun- 
damental diagram with various non-acceleration rates R 
are shown. Each point in the graph represents the result 
of 20,000 iterations starting from a single random initial 
configuration. In the diagram, we observe three distinct 
regimes of traffic flow. In between the free traffic regime 
at low density and jamming regime at high density, there 
is a distinct third regime at intermediate density. Here, 
the flux displays the sensitive dependence on the initial 
condition, which is shown in the graph as jittery struc- 
tures. 

It is certainly possible to consider the ensemble of ini- 
tial conditions among which the results could be averaged 
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FIG. 1: The fundamental diagrams of our model. The non- 
acceleration rate R is varied from R — to 0.8 as indicated 
in the graph. The maximum velocity is fixed at U — 8 for all 
cases. The size of the road is set to be L — 200. 




FIG. 2: The fundamental diagrams of two segment model 
with a common non-acceleration rates Ri = R% = are 
shown. The maximum velocity for the first segment is fixed 
at U\ =8, and it is varied as U2 =8, 3 and 1 for the second 
segment. 

out to obtain smooth curves. However, we can also ar- 
gue that jittery structure reflects a special feature of the 
model in which there can be metastable states that re- 
sult in different flux for traffics with same density but 
different initial configurations. We stress that, at this 
stage, the existence of the intermediate region is a result 
of finite size effect. This can be easily proven by the ex- 
plicit numerical calculation showing the shrinking of this 
region with larger L. At the limit L — > 00, fundamental 
diagram with two phases with discontinuous F(p) is ob- 
tained. However, this intermediate region develop into a 
distinct phase when there are more than two segments, 
S > 2 in a road. 

That fact is clearly shown in FIG. 2, in which the fun- 
damental diagram of two segment model with L = 200 
which is split into L\ = 160 and L2 — 40 is shown. Here, 
we set a common non-acceleration rate R\ = R% = and 
different maximum velocities U% and U2 for two segments. 
We assume the smaller maximum velocity for the second 
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segment, U\ > TJ%. By the inspection, we sense the exis- 
tence of analytical solutions represented by the straight 
lines connecting p = 0, p = 1 both at F = and two 
critical densities p> and p* both at the plateau flux F*. 
At the critical density p', the block of jam with Vi = U2 
and car spacing — Xi = U2 + I is formed in the second 
segment while the traffic is at full speed Vi = U\ with car 
spacing U\ + 1 in the first segment. At p*, however, the 
jam with Vi — U2 fills the entire system. With elemen- 
tary calculations, we obtain the critical densities in the 
forms 
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Li 1 

L Ui + 1 + L U 2 + 1 ' 
1 



U 2 
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We also obtain the plateau flux 



F* = 



U 2 



1 



(4) 
(5) 

(6) 



At the plateau region p^ < p < p* , the flow is limited 
by the block of jam with car spacing U2 + 1 which now 
extends into the high maximum velocity region, in which 
drivers perceive the block as bottleneck congestion. The 
length of the block is determined by the requirement that 
"free traffic" segment has the density that corresponds to 
the same flux with the block, |JBJ. 

Note that the difference between two critical densities 
p< and p* persists at the continuum limit L s — > 00 if the 
ratio L1/L2 is kept constant. This clearly shows that 
the appearance of the intermediate plateau region is not 
the artifact of the discretization nor the finiteness of the 
model. Thus this third region can be legitimately re- 
garded as the third phase. The three phase structure 
thus obtained is reminiscent to that obtained with the 
stop-and-go dynamics @- The metastable characteris- 
tics of the blocks is essential in obtaining this structure. 

The situation becomes even more interesting when dif- 
ferent non-acceleration rates are assigned to the different 
segments. In FIG. 3, the fundamental diagrams with 
a common U s and the different R s are shown. Actual 
values are chosen to be U\ = U2 = 8, R\ = 0.1 and 
U2 = 0.6. Here, in between the free traffic phase at low 
density region p < p^ and jamming phase at high density 
region p > p* , we clearly observe an intermediate phase 
where there is wild oscillations in F(p). Different out- 
come is obtained from different initial configuration. A 
most striking feature, however, is the fact that, in this 
intermediate phase, the values of the flux F is discretized 
at the fractions 



F 



V 



V+l 



V = 1,2,...,U 1 . 



(7) 



The cause of this discretization and fluctuation is again 
the appearance of metastable blocks composed of equi- 
spaced cars moving in common velocity V. The road seg- 
ment with high-i? s hinders the percolation of the blocks, 
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FIG. 3: The fundamental diagrams of two segment model 
with a common maximum velocity Ui = XJ% = 8 are shown. 
The non-acceleration rates for the first segment is fixed at 
Ri — 0.1 and that of second segment is varied as R2 = 0.4, 
0.6 and 0.8. Line connecting the data points are drawn solely 
to indicate the fluctuation between discretized flux: There 
are no points in-between the value specified by (7) in the 
intermediate density region. 



FIG. 4: The flow patterns of two segment model at p = 0.3. 
The system parameters are Ui = Ui = 8, Ri = 0.1 and R2 = 
0.5. Different initial conditions result in different metastable 
states: The graph on the left hand side shows the block of 
jam with Vi = 4, while the right hand side shows the one with 
Vi = 2. 



and helps the development of its aggregation. Thus, in 
effect, the high R s segment works very similar to the 
low-[/ s segment in previous example. A crucial differ- 
ence, however, is that the block velocity V is not preset. 
Metastable blocks with different V can appear and dom- 
inate depending on the initial condition, which results in 
non-unique value for F at each p. 

Our claim of the appearance of plural metastable bot- 
tleneck blocks are corroborated by the flow pattern shown 
in FIG. 4, which is obtained from the same simulation as 
in FIG. 3. FIG. 4(a) and FIG. 4(b) are the results of 
same parameters but different initial conditions, clearly, 
metastable congestion blocks with different internal car 
spacing (thus different block velocities) are becoming the 
limiting factors of the total flow. 

In order to further clarify the initial-condition depen- 
dence, we show the results of two different types of defi- 
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FIG. 5: The fundamental diagram of two segment model with 
a common maximum velocity Ui = U2 = 8 and different non- 
acceleration rates Ri = 0.1, R2 = 0.6. The initial condition 
is taken to be a single cluster. 
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FIG. 6: Same as FIG. 5 except that the initial condition is 
taken to be uniformly distributed, equally spaced cars. 
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FIG. 7: The fundamental diagram of two segment model 
with a common maximum velocity U\ = U2 = 8 and differ- 
ent non-acceleration rates R\ = 0.1, R2 = 0.6. The flux at 
each density is obtained by averaging over 200 different initial 
states. 



FIG. 8: The fundamental diagrams of "generic" six segment 
model with a common maximum velocity Ui = ... = Ue = 8 
are shown. The non-acceleration rates for the odd segments 
is fixed at R\ = R3 = R5 = 0.1, and that of even segment 
is varied as R2 = 0.4, R4 = 0.5, and R§ = 0.6. The lengths 
of the segments are L\ = 355, L2 = 30, L3 = 300, I/4 = 60, 
Z/5 = 110 and Lq = 45. Here, results from a single random 
initial configuration is shown for each value of the density, as 
in FIG. 3. 



nite initial states: In FIG. 5, the fundamental diagram of 
traffic flow calculated from single cluster initial states is 
shown. In FIG. 6, the initial condition is now a uniform 
distribution, meaning equally spaced cars spread over the 
whole road. In both cases, a common maximum veloc- 
ities U\ — U2 — 8 and different non-acceleration rates 
i?i =0.1 and i?2 = 0.6 are used. Even with these well- 
defined non-random initial conditions, they both show 
the fluctuating flux as the functions of density. These 
diagrams unambiguously show that the discretization is 
the results of random appearance of different metastable 
states. In FIG. 7, the result of average over different 
initial condition is shown. The system parameters are 
the same as those in FIG. 5 and FIG. 6. The flux now 
comes in between the discretized value found in the case 
of single initial condition. This is what we should ex- 
pect when the fundamental diagram is drawn from the 
accumulation of actual data from real-life traffic. 

We now move on to show that the existence of a phase 
with fluctuating discretized flux is not a specific result of 
the two-segment model, but a feature generic to multi- 
segment traffic in our model. In FIG. 8, we present the 
results of six segment traffic, which is brought up as 
a toy model of realistic single-lane "county-road" traf- 
fic. Road length is set to be L = 900. The road is 
now made up of the "normal" segments s = 1,3,5 with 
i?i = R3 = i?5 = 0.1, and the bottleneck segments 
s = 2, 4, 6, each having R 2 = 0.4, R4 = 0.5 and Rq = 0.6. 
The length of each segments are set to be L\ = 355, 
L 2 = 30, L 3 = 300, L 4 = 60, L 5 = 110 and L 6 = 45. 
These arbitrary numbers are chosen to make this case a 
generic example. A common maximum velocity U s = 10 
is used for all s. We find the result to be essentially in- 
distinguishable to the case of two segments. In a sense, 
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this is to be expected because if several jamming blocks 
are formed on a road with many segments, the one with 
smallest flux among them will eventually absorb the rest 
and becomes the limiting bottleneck, thus determine the 
total flow of the system. In this sense, analytical results 
of two segments with a common non-acceleration rate, 
(4)-(6), showing the existence of the "quantized" flux is 
the basis of all subsequent results showing the metastable 
jamming blocks. 

IV. SUMMARY AND PROSPECTS 

In real-life traffic, we often experience different mode 
of congestion on different days on same roads with sim- 
ilar number of cars. In this article, we have shown that 
this puzzling feature is captured with a very simple cel- 
lular automaton model. Our model has a property of 
self- aggregating block of jam. This jam functions as a 
dynamical bottleneck, and leads to the existence of the 
third phase at intermediate density where visually strik- 
ing non-unique flux is observed for a given density. It 
is our view, that among the multiple possible variants 
of original NS models, ones with such property are defi- 
nitely worth further attention. 

On the other hand, we want to make clear that the in- 
troduction of the alternative new rule Q is not the pre- 
requisite for the existence of the third phase. In fact, it 
can be shown through the numerical simulations that the 
three phase structure is already present in multisegment 
model with standard NS rules. With our modified rule 
that allows the formation of metastable blocks of jam, 
however, the new phase obtains distinct characteristics 
of discretized flux. 



The existence of metastable states has been already 
observed in several traffic flow models 0, 0, 01 that effec- 
tively amount to introduce some extra rules to the stan- 
dard Nagel-Schreckenberg model. In particular, early in- 
dication of discretized flux can be observed already in 
[6(. Our model is arguably the simplest of such exten- 
sions since it actually reduces the number of rules. More 
important point than the smaller number of rules is that, 
with our model, we are able to show a clear connection 
between the metastable states and the newly found in- 
termediate phase featuring the discretized flux. In light 
of our results, it is suspected that the existence of the 
metastable states observed in elaborate realistic traffic 
models do not come from detailed model-specific assump- 
tions but have a simple dynamics behind it which is cap- 
tured in our minimal model. 

The formation of metastable clusters is recently ob- 
served in a cellular automaton model of ants' trail It 
would be beneficial to compare our results with those of 
this work to gain a better insight into the dynamics of 
block of jam. 

It is worth mentioning that the state of the system 
in the newly identified third regime at intermediate 
density shares many common features with so-called 
self-organized criticality dQj. It is as though the system 
is at critical point for extended range of parameter 
values with various competing metastable states formed 
spontaneously without fine-tuning of the system. It 
should be interesting to check, for example, the exis- 
tence of 1/f type fluctuation on some measure, a hall 
mark of self-organized criticality, to further clarify the 
connection. 
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